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Vorticity generation in accretion disks around Schwarzschild and Kerr black holes is investigated 
in the context of magnetofluid dynamics derived for both General Relativity (GR), and modified 
gravity formulations. In both cases, the Kerr geometry leads to a “stronger” generation of vorticity 
than its Schwarzschild counterpart. Of the two principal sources, the relativistic drive peaks near 
the innermost stable circular orbit (isco), whereas the baroclinic drive dominates at larger distances. 
Gonsequences of this new relativistic vorticity source are discussed in several astrophysical settings. 


I. INTRODUCTION 

An exploration of the dynamics of accretion disks near 
compact astrophysical objects can advance our under¬ 
standing of phenomena as diverse as angular momentum 
transport, jet production and gamma ray bursts 
Baroclinic instability arising from misalignment between 
temperature and entropy gradients in a hot charged fluid 
rotating in an accretion disk is considered to be one of 
the most probable pathways for vorticity generation in 
astrophysics [MQ]. Such vortices (electromagnetic and 
hydrodynamic) can be amplified through several mecha¬ 
nisms such as dynamos and MRI leading to a large scale 
vortical field pTtiTS] . The resulting vortical field geome¬ 
try can be responsible for angular momentum transport 
and production, acceleration and collimation of jets in 
black holes, protostars, microquasars, etc. mm- 

In this paper, we will explore additional sources of vor¬ 
ticity generation that can contribute to angular momen¬ 
tum transport, jet production, and collimation as well 
as broaden our understanding of black hole accretion, 
in general. Unlike the traditional “baroclinic” mecha¬ 
nism, these additional drives depend on the relativis¬ 
tic effects-both special and general. Vorticity generation 
(via the relativistic drives) in the accretion disk near the 
Schwarzschild black hole was, previously, studied for a 
generalized “magnetofluid” in curved spacetime [14]. The 
relativistic drive for a pure barotropic system naturally 
emerges in the dynamics of a magnetofluid, combining 
kinematic and thermodynamical attributes of a hot fluid. 

Let us begin by recapturing the salient features of the 
magnetofluid formalism, and of earlier work on relativis¬ 
tic drives: 

l)The essence of the magnetofluid formalism (for a per¬ 
fect fluid) lies in the construction of an antisymmetric, 
hybrid tensor m 

= F^’' + {m/q)S^’', (1) 

that is a weighted sum of the electromagnetic field tensor 
(weight = charge q), and the composite (kinematic- 
statistical) fluid tensor = V^{QU^) — V^{QU^) 

(weight =mass m) [mile]. The statistical factor Q is 
the thermodynamic enthalpy. In terms of the en¬ 


tire dynamics of the relativistic hot fluid is expressible in 
the succinct equation (T is the temperature of the fluid) 




— mnV^Q 
n 


= TV^a, 


( 2 ) 


where the right-hand side is the thermodynamic force ex¬ 
pressed in terms of the fluid entropy a using the standard 
thermodynamic relation between entropy with enthalpy. 
Here, and n represent, respectively, the plasma 4- 
velocity and the number density. 

2) The 3-vector part of reduces to the more famil¬ 
iar form of 3D vortex dynamics except that the standard 
fluid vorticity is replaced by the hybrid magnetofluid vor¬ 
ticity. The addition of relativity, however, introduces 
a fundamental change; the topological helicity invariant 
Pi =< O-> (with V X ^ U being the inverse curl 
of vorticity) of an ideal nonrelativistic fluid no longer per¬ 
tains. 

Through the “distortion” of space-time, relativistic dy¬ 
namics breaks the helicity invariant even in ideal dynam¬ 
ics (cr = cr(T)); new sources and sinks appear and the 
creation and destruction of the generalized vorticity be¬ 
come possible in ideal dynamics milSlllTllIH]. Such 
sources can, therefore, be available to create vorticity in 
the accretion disk. 

3) Recently, this formalism was generalized to incorpo¬ 
rate nonminimal coupling of the magnetofluid to a curved 
background space-time m- It is quite remarkable, that 
the nonminimal coupling [introduced through fm{R), a 
function of the Ricci scalar R] changes the equation of 
motion § only minimally 

= QTV^a; (3) 


it multiplies the right-hand side with a curvature depen¬ 
dent factor Q = (1 + A/^(R)) that reduces to unity as the 
nonminimal part goes to zero, as expected. This has also 
resulted in the introduction of additional gravity-coupled 
flow field tensor which, after appropriate 3+1 decompo¬ 
sition, yields new expressions for generalized electric and 
magnetic field. The formalism, epitomized in Eqs. (§- 
Q, will be, henceforth, called Magnetofluid formalism. 

In this paper, then, we will investigate vorticity gener¬ 
ation in accretion disk for minimally as well as nonmin- 
imally coupled magnetofluid in Schwarzschild and Kerr 
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space-time. The magnitude of induced vorticity will be 
estimated in the special case when matter and space-time 
are coupled with constant Ricci scalar the simplest 
functional form of nonminimal coupling. 

We first give a summary of the derivation of gener¬ 
alized equation of motion of a new hybrid magnetofluid 
in curved background space-time. Next, the Arnowitt- 
Deser-Misner (ADM) formalism of electrodynamics [2Ql - 
[23] presented in Appendix is applied to this new for¬ 
mulation of magnetofluid and the equations obtained are 
cast into the vorticity evolution equation. These equa¬ 
tions are analyzed for accretion disks to calculate and 
estimate generalized vorticity. Finally we compare the 
relativistic drives with the more conventional baroclinic 
drive. 


II. MAGNETOFLUID FORMALISM 


Notice that the equation of motion is not yet in the 
promised “canonical” form To make progress, follow¬ 
ing Refs. nsms], we will construct the new grand unified 
vorticity tensor that reflects nonminimal coupling. 

After some patient algebra, we find that is again 

the weighted sum as , 

777 

(6) 

but with a considerably more complicated 

TD 

= (1 + xf^{R) - XRFm)S^^ + jXFmK^- (7) 

replacing — V^{QU^). We needed to 

“find” a new curvature-weighted antisymmetric flow field 
tensor 

= V^{RgU^) - \/^RgU^), ( 8 ) 


The modified theory of gravity offers an alternative 
approach to explain the inferred accelerated expansion 
of the Universe and other cosmological data by introduc¬ 
ing deviations from the GR. The so-called f{R) gravity 
modifies GR in the low energy (curvature) regime, but 
its behavior in the high energy regime has also been the 
topic of current research [241426] . To expand the scope of 
our earlier calculations, we incorporate the f{R) gravity 
[through the term Xf{R)] in the magnetofluid formalism. 
Referring the reader to Ref m for a detailed derivation, 
we simply write down, here, the two main equations: the 
modified Einstein equation (G = c = 1 ) 

-iWV’'-gnFAR) = 8nTZ,,, 

and the magnetofluid equation of motion 

(1 + A/„(i?))v^T;; = 

'qnF’' - AF„(i?)(T;; + 5 ^V)V^i?' , 

where Fg = fg{R), R is the Ricci scalar and 
is the total stress-energy tensor for both the perfect 
fluid and Maxwell’s field in curved space-time. In the 
preceding equations, r is the proper time, q is the 
charge of the particle, is the Riemann tensor and 

~ (P p)U^U^ (with = dx^^/dr) is the 

energy-momentum tensor for a perfect isotropic fluid. 
The phenomenological parameter A represents the cou¬ 
pling strength of the plasma to its background geometry, 
now modified through fm{R) and F^{R) = f^{R). The 
quantity p p = h is the enthalpy of the fluid plasma 
and often appears in the formalism as the combination 
g = hjmn with m, n, p, and p being the mass, number 
density, energy density, and pressure, respectively. 


(4) 


(5) 


rj-tflU 

total 


to derive the sought-after form. The new fluid tensor 
contains, explicitly, the coupling of flow field to gravity. 
Thus, the dynamics of a hot fluid system in curved back¬ 
ground space-time can be written into the canonical four 
dimensional vortex form 

qU,M^‘' = (1 + \fm{R))TW^cT, (9) 

the form advertised in We have “assumed” that the 
standard thermodynamic relations continue to hold; it is, 
of course, contingent upon an appropriately well-defined 
local concept of temperature in curved space-time. 

Notice that, when A = 0, reduces to its minimally 
coupled counterpart tensor defined in [HUS]]- 

Equation is the main result that describes the mag¬ 
neto fluid dynamics in curved space-time. It reveals 
that a charged relativistic fluid, coupled nonminimally 
to gravity, obeys a 4D vortex dynamics like its grav¬ 
ity free and minimally coupled (to gravity) counterparts. 
The new grand vorticity tensor subsumes earlier limiting 
cases in a transparent manner. 

We would now apply the above formulation to inves¬ 
tigate the vorticity generation in accretion disks around 
black holes. To do calculations in terms of familiar quan¬ 
tities, an appropriate 3+1 decomposition of the space- 
time is necessary; it is presented in Appendix Next, 
we present the 3D vortical dynamics in order to facilitate 
computation of vorticity generation. 


A. Vortical dynamics 

With the 3+1 decomposition presented in Appendix[Al 
the spacelike projection, i.e., 7 ^ ^ projection of the uni¬ 
fied field equation of motion gives us the momentum 
evolution equation 


aqV^ + qr{v x U) = -(1 + Xfm{R))TVa 


( 10 ) 
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whereas the timelike ( n^) projection gives the equation 
of energy conservation 


aqVv-^ = T(1 + Xfm{R)){Cta - (3 ■ Vct), 


( 11 ) 


where ^ and O are, respectively, t he ge nerali zed electric 
and magnetic fields given by Eqs. (A 6 ) and (A7) in Ap ¬ 
pend ix [A} Also, Q^, r, and v are defined in Eq. 1 ®,®, 
and (A3). Sources responsible for magnetic field gener¬ 
ation, in particular, the sources that are gravity driven, 
can be derived from the generalized vorticity evolution 
equation (which is really the generalized Earaday’s law) 
by manipulating Eq. ( [1Q| ). 

Since is an antisymmetric tensor, the divergence 
of its dual is zero, i.e., = 0. Taking the 7 ^ ^ 

projection of the preceding identity, we derive 


Ct^ = - V X {ai) - aOn, 


( 12 ) 


where C denotes Lie derivatives with jCt = dt along 

= [/3,fl], and the expansion factor 0 is defined in 
Appendix 

It should be noted that, even in the absence of nonmin- 
imal coupling to gravity (A = 0), (minimal) coupling to 
gravity still manifests in the formalism. Equation (12), 


in conjunction with Eq. (10), gives us the vorticity evo¬ 
lution equation 


£tfi-Vx(i;xj1)-7i2+a0j1 = Vx (^T(i + A/„(i?))Vcr 

(13) ' 

All terms on the left-hand side operate on the vortic¬ 
ity 3-vector Q while the right-hand side provides, just as 
in the conventional picture, possible sources for vorticity 
generation. The left-hand side, however, has lot more 
structure than the conventional 3D vortex dynamics; the 
first two terms reflect the standard Helmholtz vortical 
dynamics, while aGfl and are nontrivial gravity 

modifications. Thus, the gravity coupling does, funda¬ 
mentally, modify the projected 3D vortex dynamics, in 
spite of the fact that the 4D vortex equations had exactly 
the same form. 


III. VORTICITY GENERATION 

To apply the formalism to vorticity generation in as¬ 
trophysics, specifically in accretion disks around compact 
objects like Schwarzschild and Kerr black holes, we have 
to specify the space-time geometry the space-time metric 
that controls the motion of plasma particles. The stan¬ 
dard metric describing the stationary and axially sym¬ 
metric (or spherically symmetric) spacetime for Kerr (or 
Schwarzschild) black holes can be written as [27] 

ds'^ = gudf + 2gt4,dtd(t)+grrdr‘^+9e9d0‘^+g4,4,d(t)^. (14) 

The exploration of geodesic motions of plasma in ac¬ 
cretion disks will allow us to compute various relevant 


physical quantities. Since we are interested only in the 
timelike geodesics in thin accretion disks, the Euler- 
Lagrange equations can be derived from the Lagrangian 
for the above stationary and axisymmetric spacetime, 
2L = ds^ jdr^ = —1, with r being the proper time 
along timelike geodesics. Thus, the corresponding Euler- 
Lagrangian equations describing the timelike geodesics in 
the equatorial plane take the form m) 

dt _ + Lgt(f) 

7 2 ’ 

dr gf^ - gu94>4> 


d(j) 

dr 


Egtcj) + Lgtt 

~~2 ’ 

dtcf) ~ 9tt9(f)(f) 


(16) 


(dr\ E‘^9cf)cf) + ‘^ELgtcf) + L‘^gtt .. 

drr (zr I = -^^- = ^e//, 

where E and L are specific energy and specific angu¬ 
lar momentum respectively. Eor stable circular orbits in 
the equatorial plane, using Veff = 0 and dVeff/dr = 0, 
the constants of motion including angular velocity uj are 
found to be 


^ _ gtt P gtcf)^ 

(18) 

'sj gtt ‘^gtcj)^ g(j)(j)^‘^ 

P _ gtcj) + g(f)cf)^ 

(19) 

\l gtt ‘^gtcp^ gcf)cf)^‘^ 

gtcf),r T \/ {gtcf),r)‘^ gtt^rgcfxp^r 

( 20 ) 

g4>4>,r 


The Lorentz factor for particles can be derived from 


P_ ^gcpcf) T Egf(i) _ 1 (21) 

9 tcj) gttgcjxcj) gtt ‘^gtcf)^ ^ 00 ^^ 

Moreover, for any given scalar function P, the gradient 
is defined as 

VP = —9^P H— -^^OqP cq -\— - d(f)P e^. ( 22 ) 

x/gEr xj goo x/gE^ 

Next, we assume a thin accretion disk with zero lati¬ 
tudinal speed vq = 0 for the plasma; we will also assume 
that the radial velocity of the plasma is negligible com¬ 
pared to the orbital velocity V(f) » Vr. The orbits of the 
plasma constituents are also taken to be almost circular 
(r ^ 0 ). Since our formalism is based on perfect fluid, we 
can also assume the plasma to be barotropic with its pres¬ 
sure depending on density only, i.e., a = E{T). Then, 
with this assumption, we can write the relation between 
the temperature and entropy gradient as TVcr = y/c^VT, 
which evidently will cause the baroclinic drive Nb to van¬ 
ish, where x is a dimensionless quantity of order unity. 

To compute the appropriate temperature profile in the 
region of interest in the accretion disk, we follow the pre¬ 
scription presented by Novikov and Thorne [28]. It turns 
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out that, for M = 14.3M0 and M = 0.472 x 10^^gs~^^ 
the inner- and outermost stable circular orbits (average 
width of the accretion disk) are located mostly in the op¬ 
tically thick region of the accretion disk ; this is true for 
both geometries. Therefore, for an optically thick region, 
using the Stefan-Boltzmann law, the temperature profile 
can be written as [26l |27] 

T(r) = ^ , (23) 

\0'SB J 

where asB is the Stefan-Boltzmann constant and 2 ; is 
redshift due to gravitational effects. For 0 = 7r/2, and 
for a vanishing disk inclination angle, the redshift can be 
written as 1 + 2 : = F. Here f(r) is the energy flux for a 
relativistic accretion disk presented in Refs. [[26l[29]] by 
Page and Thorne as 



f(r) = 


i\To uj ,7. 

{E - ujiy 



(24) 


FIG. 1: Temperature profile for the Schwarzschild (blue) aud 
Kerr black hole (red) from x = 1.5 to x = 30. 


where Visco is the radius of the innermost stable circular 
orbit in the accretion disk, and Mq is the mass accretion 
rate. The temperature and the Lorentz factor profiles, 
displayed, respectively, in Figs 0 and Q, reveal simi¬ 
lar general features for the Kerr space-time: increasing 
from their corresponding value at r^sco? they reach a peak 
at some radius and then monotonically decrease as we 
move away from the center of the corresponding black 
holes m- However, only the temperature profile in the 
Schwarzschild geometry shares the similar feature. These 
features of the temperature and the gamma (F) profiles 
will manifest in the vorticity generation as well as in the 
relative strength between the corresponding relativistic 
and classical drives. 

Throughout this paper, we will use parameters ob¬ 
tained from observation on the Galactic black hole 
Cygnus-XI as representative: M = 14.38M0, M = 
0.472 X and a= 0.99, where M and a are, re¬ 

spectively, the mass, and angular momentum per unit 
mass of the black hole [SilISS- Figure!^ shows the tem¬ 
perature profile in Kerr (Schwarzschild) geometry from 
r = l.bvg (r = Gvg) to r = 30r^ with Vg = GMjc^ in 
the normalized unit of x = r/vg. The profile shows a 
peak temperature between 10^ — 10^K which drops as 
we move from the event horizon radially outward. These 
profiles are used in this paper to calculate the vorticity 
generation in the accretion disk. 

A. Schwarzschild geometry 



FIG. 2: Lorentz factor for Schwarzschild (green) and Kerr 
black hole (red) from x = 1.5 to x = 30. The dashed 
line (blue) indicates inner-most circular orbit (isco) for the 
Schwarzschild blackhole. 

the time translation symmetry leading to a vanishing ex¬ 
trinsic curvature, the new term involving 0 on the left 
hand side disappears. Thus, the structure is precisely 
like the 3D vortex dynamics. The simplified vortical evo¬ 
lution equation can be used to approximately compute 
the weak field seed generation in the hot fluid system in 
the accretion disk in Schwarzschild geometry. 

The relevant space-time metric elements are 


For a spherically symmetric and static space-time 
(Schwarzschild space-time), the above vortical evolution 
equation ( p^ reduces to the one presented in Ref. El, 

i.e., — Vx(Fxf1) = Vx . Since the 

spherically symmetric and static space-time can be foli¬ 
ated without any shift function /3, and the foliation obeys 


2 /-I ^0 X 

gtt = -a = -(1-) ; 


gee = r ; 


= (i-7rt 

i = sin^ 0. (25) 


Then, using the Eqs. ( TsjpT ), we can calculate the 
orbital velocity and gamma factor for the orbital motion 
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of plasma in the accretion disk. Since the radial velocity 
of the plasma is assumed to be negligible compared to 
the orbital velocity V(fy » the orbits of the plasma 
elements are taken to be almost circular (r « 0). Then, 
inserting the Schwarzschild metric elements, ( [^ , in Eqs. 
( T5p7| ) and imposing Veff = 0 and dVeff/dr = 0 reveals 
that there exists one stable circular orbit at x > 6 and 
one unstable circular orbit at 6 > x > 3 [[22]]. This 
also dictates the applicability of the temperature profile 
in accretion disk. The temperature profile used in this 
paper is valid in the region of stable circular orbits. 


B. Kerr geometry 


For an axisymmetric (but not spherically symmetric) 
stationary system, like the Kerr black hole, our previous 
assumption of the zero shift function, = 0, is no longer 
valid. Consequently, the pertinent equation (13), in 
general, does not show any similarity to standard 3D 
vortex dynamics. The shift function for rotating black 


holes can be taken to he (3 = —uj with respect to a 
zero angular momentum observer. The term involving 
the shift function, however, will give zero contribution 
since we assume that both uj and ft have only radial de¬ 
pendence. In addition, it can be further shown that the 
term involving the expansion factor 0 vanishes. Thus, 
the vorticity evolution equation, even for Kerr geometry, 
will resemble the standard 3D vortex dynamics. 


The relevant space-time metric elements in Boyer- 
Lindquist coordinates are [26] 


„ _ „ 2 _ (.^r - _ 

9tt ^ ^2 ’ drr 

gt4, = +a^ - Ar) ; = 


A/ 

(26) 


where = (r‘^ o?) — 2r^r. Then, again inserting the 

Kerr metric elements, ( [^ , in Eqs. ( T5][T7 ) and imposing 
Veff = 0 and dVeff/dr = 0 reveals that the inner most 
circular orbit in the Kerr black hole is located at x = 
1.4545 for a = 0.99r^, which was taken into account in 
deriving the temperature profile for the accretion disk in 
Kerr geometry. Note that we will assume a = 0.99r^ 
throughout the rest of the paper. 


C. Computing vorticity 

For both Schwarzschild and Kerr configurations, com¬ 
putation of vorticity generation requires knowledge of the 
(p dependence of the temperature profile. However, the 
most commonly used temperature profiles (including the 
GR corrected ones) for accretion disks show only radial 
dependence. Previously, an estimate of vorticity genera¬ 
tion was computed using an average temperature of the 


accretion disk [14]. However, as shown in Fig. 0 , Gen¬ 
eral Relativity restricts the application of an average disk 
temperature throughout the disk as it involves regions of 
unstable orbits leading to nonlinear behavior. Moreover, 
plasmas orbiting in accretion disks for both black hole 
configurations undergo gravitational radiation reaction, 
which for a Kerr black hole can cause a plasma particle 
to lose as much as 42 percent of its initial energy as it 
approaches the event horizon [22]. 

Toroidal temperature dependence is created due to the 
gravitational radiation by the orbiting plasma particles; 
the induced radiation reaction, in turn, makes the sta¬ 
ble circular orbits deviate slightly from geodesic motion 
[22]. A particle, initially in a circular orbit at x > 6 
{x > 1.4545) for the Schwarzschild (Kerr) metric, slowly 
spirals into smaller nearly circular orbits as it radiates en¬ 
ergy until it reaches the orbital radius x = 3 {x = 1.4545), 
where the orbit becomes unstable. Therefore, the stable 
circular orbits do not close in either geometry, and, de¬ 
pending on the magnitude of radiation reaction, the spa¬ 
tial orbital trajectory in the equatorial plane is assumed 
to be represented by r = r(0), a solution to the geodesic 
equation relating coordinates r to (j). 

Next, we assume, without loss of generality, that the 
orbits of the spiraling plasma elements in the disk can be 
approximated by r(0) = where vq is the initial 

radial distance of the plasma particles from the center of 
the black hole, and ( is the parameter that controls how 
tightly a nearly circular orbit spirals around the black 
hole. In general, the factor ( can be a complicated func¬ 
tion of black hole mass as well as the energy and the an¬ 
gular momentum of plasma elements. Determined from 
the geodesic equation relating coordinates r to 0 with ap¬ 
propriate boundary conditions, can be a function of the 
radial distance. However, since our focus is on the region 
of accretion disk over which the timelike geodesic orbits 
of the plasma elements are closely bound, we can assume 
the spiraling parameter C to be a constant. A rapidly 
varying spiral with varying C will contribute more to the 
vorticity as can be seen from Eq.(27). 


As mentioned earlier, the radiation reaction, can be the 
source for slight deviations from the closed stable circular 
orbit, thereby imparting temperature variations around a 
spirally circular orbit. Thus, for the spiral orbit, r(0) = 
dT/d(\) = (dTjdr)(drjd(\)) = —C,r(dTjdr), and 
the seed vortical field \(}\ may be estimated as 


fi(r) = - 


xhcrC 


1 


0. sj 9rr9cf)(f) 


dr{T~^){drT)M (9, (27) 


where At is the characteristic time for linear vorticity 
generation under which the changes in space-time geom¬ 
etry are negligible. Therefore, we choose this time-scale 
to be At = 271 and, as expected, this coordinate time 
interval At is related to the proper time interval Ar by 
At = TAr. For an observer far away from the accre¬ 
tion disk under observation, these two time intervals are 
practically the same. For the Schwarzschild geometry. 
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FIG. 3: Vorticity magnitude (in the units of Gauss) for the 
Schwarzschild black hole with different ( values between 6.0r^ 
and 30r^ for M = 14.SM© 



FIG. 4: Vorticity magnitude (in the units of Gauss)for the 
Kerr black hole with different ( values between 1.5r^ and 30r^ 
for M = 14.SM© with an accretion rate M = 0.472 x 10^^^/s 
and a=0.99. 


Eq. (27) simplifies to 


following metric components 


nix) = §, (28) 

qV^Vg 

X = r/rg is the normalized distance, and g = —e is the 
electron charge. 

Figure and 1^ show the radial dependence of the mag¬ 
nitude of the vorticity generated in the accretion disk 
plasmas embedded in Schwarzschild and Kerr geometries 
for three different choices of Both figures show a dras¬ 
tic reduction in \ fl\ as the temperature maxima are ap¬ 
proached (temperature gradient going to zero), and then 
|0| picks up as we go over to the other side of the max¬ 
imum. In addition, the existence of a second dip in the 
Kerr vorticity profile can be attributed to the vanish¬ 
ing of the gradient of the corresponding Lorentz factor. 
Therefore, the effect of the relevant gamma and tem¬ 
perature gradients indicates a distinct vorticity profile in 
Kerr space-time, where vorticity changes direction twice 
before it gradually decays radially outward in the disk. 

However, as we move radially outward in the disk, 
both induced vorticities decrease almost at the same rate 
maintaining their difference. 


D. Black hole accretion disk in modified gravity 


The results of last section were derived for the ac¬ 
cretion disk fluid minimally coupled to standard GR. 
We will, now explore the changes in vorticity generation 
brought about by modified gravity, more precisely, by the 
simplest functional form of f{R) = Ro = constant. Fol¬ 
lowing the discussion presented in Sec. |III[ plasma parti¬ 
cles in the Schwarzschild accretion disk are described by 


9tt 



Qrr — 9 5 


= siM 6>, 

(29) 


where Rq is the constant Ricci scalar. First, setting 
= 0? we obtain the condition on Ricci scalar Rq 

Uox^ - 12x + 24 = 0 (30) 


with Rq = RoVg and x = rjcg. To get a black hole with¬ 
out any naked singularity and with event/cosmological 
horizons, the solution of Eq.(30) restricts the value for 
IZo to be (—00,4/9) ([H])- Eurther restrictions on Rq 
to isolate the geodesics for stable circular orbits can be 
imposed by demanding W// = 0, dVeff/dr = 0, and 
£Veff/d^r>0, 


12(6-Xe) 

(15 — 4xc)x3 


(31) 


with Xc being the radius of the stable circular orbit. The 
above relation also reveals that there exist one innermost 
and one outermost stable circular orbit, and the upper 
limit on Rq reduces to 2.85 x 10“^. However, the choice 
of an exact value of Rq befitting our current analysis 
depends on the average width, the temperature, and lu¬ 
minosity profiles of black holes. Thus, if the temperature 
and luminosity profiles prescribed by Page and Thorne 
[29] are taken into account, it turns out that, for f{R) 
Schwarzschild black holes, the new range of Rq further 
reduces to (—oo,10“^]. Then the only judicious choice 
turns out to be Rq = 10“^, for it satisfies the average 
radius of the outer edge of a Schwarzschild black hole 
accretion disk at r ~ 70r^ by setting the innermost and 
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outermost stable circular orbits, according to Eq. (31), 
at Xc = 6 and Xc = 143.45 [26] . 

A similar analysis, carried out for the Kerr metric in 
modified gravity with the metric elements. 


2 {Ar - a^) 

9tt = -a = -— 


9t4> = - ^r) ; 

~ "2^2 


A/ 


(32) 


with 


^0 2 \ 
12 ^ 


- 2r^r, 

(33) 

Rq 2 

12 ^ ’ 

(34) 


demands IZq (upon setting 1 / g^r = 0) to satisfy 




1Zqx‘^ 

12 


2x = 0. 


(35) 


While Eq.(35) yields the range for IZq G [—0.3,0.6], a 
new range for TIq G (0, 0.6] emerges if a Kerr black hole 
with two event horizons and one cosmological horizon 
is demanded. However, once again, upon demanding 
Veff = 0, dVefffdr = 0, and ^V^ff/^r > 0 along 


with the appropriate temperature and luminosity pro¬ 
files mentioned above, we find the stable circular orbits 
can exist only for T^o ^ [—1.2 x 10“^, 6.67 x 10“^]. To 
maintain the consistency in our numerical plots, we again 
choose 7^0 = 10“^, for it satisfies the the average radius 
of the outer edge of a Kerr black hole accretion disk at 
r ~ 16r^ by setting the innermost and outermost stable 
circular orbits at x^ = 1.4545 and x^ = 143.45, respec¬ 
tively. It should be noted here that temperature profiles 
used for both classes of accretion disks remain the same 
as long as we choose T^o ~ [2S]- 

The general expression of vorticity generated in modi¬ 
fied gravity is 


fi(r) = -(1 + [Q^T) At 0, 

Q. y 9rr9(f)(f) 

(36) 


with Vm and Xfm{R) = denoting the modified 
Lorentz factor and the nonminimal coupling of plasma 
to f{R) gravity, respectively. Similarly to Eq.(28), we 
have an analytical expression for modified Schwarzschild 
spacetime 


= + (37) 

qy^XTg 

where am is associated with modified Schwarzschild met¬ 
ric component gu- 

Eigures and show the plot of this generalized vor¬ 
ticity \Q\ (in the disk) as a function of the distance 
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FIG. 5: Vorticity Q (in the units of Gauss) for Schwarzschild 
black hole with different ^ values in modified gravity within 
6.0r^ and 30r^ for M = 14.SM©. 



FIG. 6: Vorticity Q (in the units of Gauss) for Kerr black hole 
with different ( values in modified gravity within 1.5r^ and 
30r^ for M = 14.SM© with an accretion rate M = 0.472 x 
10^^g/s and a=0.99. 


X = r/rg in the gravitational field of a black hole of mass 
M = 14.SM© and TZq = 10“^. The kinklike behavior of 
11 vs X = r/vg, as earlier, originates in the vanishing 
gradients of temperature and the Lorentz factor. Differ¬ 
ent values of C, used in the Figs. <§ and capture 
the influence of the modified gravity on toroidal temper¬ 
ature fluctuations. Here, it should be emphasized that 
vorticity magnitude profile in the modified gravity can 
be drastically different for arbitrary functions of /(i?), 
which signifies how the matter is coupled to its back¬ 
ground spacetime. 

It must be noted that, as maximum temperature ob¬ 
tained in f{R) models turns out to be lower than Cygnus 
X-1 observation values, the probability of the existence 
of an accretion disk in the f{R) Schwarzschild black hole 
becomes very slim m- Still we present the result here for 






























f{R) Schwarzschild black hole as an analytical example. 

E. Relative strength of relativistic and baroclinic 
drive 

The generalized magnetofluid dynamics, derived for 
a barotropic equation of state, will have only the rela¬ 
tivistic vorticity drives. We would now like to compare 
the relative magnitudes of the relativistic, and possibly 
baroclinic, vorticity sources. At first glance, it is evident 
(from the metric) that the relativistic source will domi¬ 
nate the baroclinic source as we approach the event hori¬ 
zon. But to make a comparison between the relativistic 
and baroclinic drives, let us construct a simple baroclinic 
drive of the form VT x Va where 

e is a measure of the departure from strict barotropic 
behavior. The preceding approximation introduces a 
smaller toroidal temperature variation that quantifies the 
nonbarotropic component. 

In Figs. 0 and 0, we compare the relative strength 
of the relativistic and baroclinic drive for black hole ac¬ 
cretion disks in Schwarzschild and Kerr geometries. We 
plot the relative magnitude as a function of the radial 
distance, x = r/r^; both plots start from their respec¬ 
tive isco. In both cases, the relativistic drive becomes 
dominant as we approach the inner most stable orbit. 
For smaller values of e, we see the magnitude of the 
relativistic drive keeps increasing as the departure from 
barotropic fluid is minimal. Also, the dashed line shows 
the ratio to be unity and from both figures we see that 
the relativistic drive remains dominant, for smaller val¬ 
ues of e; the dominance continues to longer distances 
from innermost stable circular orbit. In both cases, the 
Lorentz factor plays an important role in determining 
their relative strength. Thus, the sudden dip in the rela¬ 
tive strength profile in Kerr metric can be attributed to 
the vanishing gradient of the corresponding Lorentz fac¬ 
tor inherent in the relativistic drive. Moreover, the rela¬ 
tive magnitude in the Kerr black hole is less than that in 
Schwarzschild black hole because of the significant energy 
loss of plasmas in Kerr space-time due to gravitational 
radiation reaction as it approaches the nonstable orbits. 
The jump in relative intensity from e = 0.01 to e = 0.05 
remains significant in both geometries. 

IV. CONCLUSION 

We have explored in this paper the possibility of gen¬ 
erating what may be called “generalized vorticity” in ac¬ 
cretion disks surrounding compact gravitating objects, 
in particular, the Kerr and Schwarzschild black holes. 
The accretion disc plasma is coupled nonminimally [via 
f{R) gravity] to the surrounding specified space-time. 
Although vorticity can be generated by the well-known 
baroclinic mechanism ( nonzero VT x Vcr), we have con¬ 
centrated, in this paper, on exploring what are classed as 



FIG. 7: Ratio between relativistic and baroclinic drive in 
Schwarzschild black hole with e=0.01 (black), 0.05(green) and 
0.1 (blue).The red line represents the ratio of unity. 



FIG. 8: Ratio between the relativistic and baroclinic drive in 
the Kerr black hole with with e=0.01 (black), 0.05(green) and 
0.1 (blue). The red line represents the ratio of unity. 


relativistic drives stemming from space-time distortions 
caused by special as well as general relativistic effects. 
We find that the Kerr geometry, due to the intrinsic rota¬ 
tion, is a more efficient vorticity generator as compared to 
the Schwarzschild counterpart. We also observe a slight 
increase in vorticity generation even if an extremely weak 
f{R) (nonminimal) coupling, f{R) = Rq = constant^ 
is turned on, which implies that a noticeable change in 
vorticity generation will be observed in strong f{R) cou¬ 
pling. Physically, the increase in efficiency is directly re¬ 
lated to the amount of deviation from the circular orbits 
(of the plasma particles) caused by the distorted geome¬ 
try. 

We also compared the efficiency of the relativistic drive 
with a model baroclinic drive. The strength of the baro- 
clinic drive is given in terms of a parameter e, that mea- 
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sures the departure from strict barotropy. For most ac¬ 
cretion disk plasmas in quasiequilibrium, e is expected to 
be small. We find that, for reasonable values of e, the 
relativistic drive is dominant for the disk regions nearer 
to the compact object; as the distance increases, and for 
relatively larger e, the barcolinic term becomes compara¬ 
ble to the relativistic drive. This vorticity growth occurs 
in orbital time scale, similar to MRI instability, which 
later can amplify through several dynamo mechanisms 
[M35]. 

Unlike GR-MHD, our formalism can be extended to 
multifluid species with each species obeying its its own 
vortical dynamics j36l [37] . In addition, the induction 
equation in this formalism has source terms for vortic¬ 
ity generation; these sources can catapult the accretion 
disk to a state of finite vorticity from one with no vor¬ 
ticity (electromagnetic and hydrodynamic). This formal¬ 
ism can also be studied in the context of vortex gener¬ 
ation in the protoplanetary disks near a supermassive 
star where gravity plays a dominant role. Hyrdody- 
namic simulations indicate these proto-planetary disks 
to be inherently baroclinic due to negative radial en¬ 
tropy gradient [71 [131 ES] • The turbulence caused by the 
baroclinic instabilities is found to be suitable for angular 
momentum transport and vortex formation in the disk 
which are suggested to lead to planet formation. Our 
model, if applied to the evolution of the protoplanetary 
disk, will provide additional sources for angular momen¬ 
tum transport and vortex formation even in a barotropic 
disk. 

Regardless of the details, the curvature drive (in mini¬ 
mal as well as in nonminimal gravity) will always gener¬ 
ate “generalized” vorticity which, either by itself or after 
amplification through a generalized dynamo mechanism, 
can provide a mechanism for angular momentum trans¬ 
port. These amplified vortical fields can also collimate 
the jets emanating from the disk where field lines coro¬ 
tate with the disk, by the flux freezing theorem. The 
plasma leaving the disk can drag the field lines due to its 
large conduction coefficient which wraps the field lines 
around the rotation axis. The field lines then exert a ra¬ 
dial force which can compress the jet of plasma leading 
to jet collimation [2]. 
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Appendix A: 3-|-l Dynamics of GravitoMagnetofluid 

The approach chosen for the 3+1 splitting selects a 
family of foliated fiducial 3-dimensional hypersurfaces 
(slices of simultaneity) labelled by a parameter t = 


constant in terms of a time function on the manifold. 
Furthermore, we let t^ be a timeline vector whose in¬ 
tegral curves intersect each leaf of the foliation pre¬ 
cisely once and which is normalized such that t^V^t = 1. 
This t^ is the ‘evolution vector field’ along the orbits of 
which different points on all E^ = E can be identified. 
This allows us to write all space-time fields in terms of 
t-dependent components defined on the spatial manifold 
E^. Lie derivatives of space-time field along are iden¬ 
tified with “time derivatives” of the spatial fields since 
Lie derivatives reduce to partial time derivative for an 
adapted coordinate system t^ = (1,0, 0,0). 

Moreover, since we are using the Lorentzian signature, 
the vector field t^ is required to be future directed. 
Let us decompose t^ into normal and tangential parts 
with respect to E^ by defining the lapse function a and 
the shift vector /3^ as t^ = an^ + with p^n^ = 0, 
where n^ is the future directed unit normal vector 
field to the hypersurfaces E^. More precisely, the 
natural timelike covector n^ = (—a, 0,0,0) = —aV^t 
is defined to obtain n^ = {1/a^ /a) which satisfy 

the normalization condition n^n^ = —1. Then, the 
space-time metric induces a spatial metric 

by the formula + n^n^y. Finally, the 3+1 

decomposition is usually carried out with the projection 
operator q^ ^y = ^ n^n^y, which satisfies the 

condition = 0. Also, the acceleration is defined as 

Now, with the above foliation of space-time, the space- 
time metric takes the following canonical form m 

ds^ = —a^dt^ + ^ij{dx^ + P^dt){dx^ + p^dt)^ (Al) 

and it immediately follows that, with respect to an Eu- 
lerian observer, the Lorentz factor turns out to be 

T=[a^- jij + 2/3^+ v^v^)] , (A2) 

satisfying dr = dt/T, where u* is the ith component of 
fluid velocity v — dxjdt. Then the decomposition for the 
4-velocity is m 

U>^ = aTn>^ + Tj'^ ^v’', (A3) 

with n^U^ = —aV. 

Now, since our unified anti-symmetric field tensor 
is constructed from the antisymmetric tensors and 
we apply the ADM formalism of electrodynamics 
presented in [2QH23] to define the generalized electric and 
magnetic field, respectively, as 

r = (A4) 

and thus express the unified field tensor 

(A5) 

We remind the reader that the generalized magnetic field 
and the generalized vorticity are essentially synonymous. 
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Using the definition of the unified field tensor the 

expressions of 3D generalized electric and magnetic fields 
turn out to be 




A/„(i?) - XRF^{R))V{agT) 


- -XFrr,iR)ViagRr) 


TY) 

-(l + A/„(i?)) 

Tfl 

— {l + Xfm{R) 

qa 


2a • (G^v) + -QG^v 
o 


- —AF„ 

qa 


XRF^iR)) (^CtiGTv) - C^igVv 
(R) (^CtiGRTv) - C^igRTv 


(A6) 


X = B 


m 


(1 + \fm{R) - \RFm{R))V X {GVv) 


^\Fm{R)-V x{RGTv), (A7) 


where a = a^ and 0 are, respectively, the shear 
and expansion of the congruence, defined as = 

~ ® We have also used 

the relation V^Uy = —ayn^^a(^i3-\-^0^^y to derive (A6). 
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